Slow decay of magnetic fields in open Friedmann 

universes 



John D. Barrow 1 and Christos G. Tsagas 2 
1 DAMTP, Centre for Mathematical Sciences 
University of Cambridge, Wilberforce Road, Cambridge CB3 OWA, UK 
2 Section of Astrophysics, Astronomy and Mechanics, Department of Physics 
Aristotle University of Thessaloniki, Thessaloniki 54124, Greece 



Abstract 

We describe how magnetic fields in Friedmann universes can experience superadiabatic 
growth without departing from conventional electromagnetism. The reason is the relativis- 
tic coupling between vector fields and spacetime geometry, which slows down the decay of 
large-scale magnetic fields in open universes, compared to that seen in perfectly flat models. 
The result is a large relative gain in magnetic strength during the pre-galactic era that can 
lead to astrophysically interesting _B-fields, even if our universe is only marginally open today. 

PACS numbers: 98.80.-k; 98.62.En; 98.65.Dx 

Magnetic fields in Friedmann universes are widely believed to decay adiabatically regardless 
of the electrical properties of the cosmic medium. Consequently, large-scale £?-fields are expected 
to dilute as a~ 2 , where a is the cosmological expansion scale factor. This widespread perception 
has its roots in the conformal invariance of Maxwell's equations and the conformal flatness of 
the Friedmann-Robertson- Walker (FRW) spacetimes. The two are thought to guarantee that 
the rescaled magnetic vector B a = a 2 B a evolves as in Minkowski space. This then ensures that 
the magnetic flux remains conserved and consequently that B a cx a~ 2 irrespective of the electric 
properties of the universe. Following [I], conformably flat spacetimes, like de Sitter space and 
the FRW models, can be written as time-dependent rescalings of Minkowski space. Then, the 
conformal triviality of Maxwell's theory guarantees that, when written on an FRW or a de 
Sitter background, the wave equation for the n-th Fourier mode of the rescaled magnetic vector 
B a = a 2 B a takes the Minkowski- like foruQ 

£' ( ' n) +n 2 £ (n) =0, (1) 

with the primes indicating conformal-time derivatives. The above guarantees that B a oc a~ 2 and 
therefore an adiabatic decay for the l?-field. Hence the belief that to modify the a~ 2 -law and 
achieve a superadiabatic-type of magnetic amplification we need to abandon either the FRW 



1 To avoid confusion, we note that B a is the magnetic vector, B a = a 2 B a is the rescaled magnetic vector, B( n ) 
and B/ n ) are their associated Fourier modes and pB = B 2 /8iv (with B 2 = B a B a ) is the magnetic energy density. 
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models or conventional electromagnetismo The usual choice is to follow the latter route [3], and 
the literature contains a plethora of mechanisms that slow the adiabatic depletion of magnetic 
fields down by departing from standard electromagnetic theory (see [1] for a representative 
though incomplete list). This is not always necessary, however, because the argument leading 
to Eq. ([JJ - and the equation itself - holds only when the FRW model has flat spatial sections 
(or when the space is exactly de Sitter). 

All three Friedmann universes are conformally flat but they are not identical. Differences in 
the geometry of their 3-spaces ensure that the conformal factor of the spatially curved models 
has an additional spatial dependence and therefore it no longer coincides with the cosmological 
scale factor. The associated line elements have the general form [51 [6] 



where a = a(£,R) is the conformal factor and dQ 2 = d6 2 + sin 2 6d(j) 2 (see [6] for details - par- 
ticularly on the open-FRW case). The above shows that Friedmann models with non-Euclidean 
spatial geometry cannot be written as simple, time-dependent rescalings of Minkowski space. 
For our purposes, this is the key difference between the flat and the rest of the FRW cosmolo- 
gies. Putting it in geometrical terms, there is no global one-to-one correspondence between 
curved FRW models and Minkowski space: the conformal transformations mapping the associ- 
ated spacetimes are only local [5] . Consequently, a Minkowski-like evolution for B a is not a priori 
guaranteed in these models. Rescaling the magnetic field with a space-independent conformal 
factor does not work in the case of non-zero 3-curvature, as it does on a spatially flat FRW 
background, and the associated wave equation need not take the form of JT]). Thus, the rescaled 
/3-field should show a Minkowski-like behaviour only locally (i.e. on small scales). On large 
scales, where the curvature effects are important, one would in principle expect to see deviations 
from the standard B a oc a _2 -law. 

This is indeed what happens. Consider, for example, the simple case of a source-free elec- 
tromagnetic field on a general FRW background. Inflation is believed to generate such clas- 
sical electromagnetic fields, by stretching the associated small-scale quantum fluctuations to 
super-horizon scales. In the absence of sources, the magnetic component of the Maxwell field 
propagates according to the linear wave equation [7] 



Here, B a is the magnetic vector measured in a frame moving with 4- velocity u a . The latter 
defines the comoving (fundamental) observers and it is normalised so that u a u a = —1. Also, 
the scalars p and p represent the matter density and pressure respectively, while lZ a b is the 
background 3-Ricci tensor (see [8] for further details). Finally, overdots indicate proper-time 
derivatives, D 2 = h ab V a Vb is the 3-D Laplacian (with h a b = g a b + u a Ub projecting orthogonal to 
u a ) and H = a/a is the Hubble parameter. Expression ([3]) can be obtained either by linearising 
the nonlinear equation (40) given in [7] (see also Eq. (45) in the same paper), or by simply 
recasting the first-order relation (28) of [9]. In addition, the electromagnetic field vanishes in 

2 Superadiabatic amplification is a concept originally introduced in gravitational wave studies [2] and usually 
implies a magnetic decay-rate slower than the adiabatic one and not actual amplification. 



ds 2 = a 2 (£, R) [-df + dR 2 + R 2 dn 2 ] 





(3) 
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the unperturbed FRW background, which frees our analysis from gauge-related ambiguities. The 
key quantity for our purposes in Eq. ([3]) is the magneto-curvature term, lZ a bB b , which results 
from the non-commutativity of covariant derivatives in non-Euclidean spaces. It reflects the 
fact that vector sources, like the Maxwell field, 'feel' the curvature of space through the Ricci 
identities - in addition to Einstein's equations. 

Introducing the rescaled B a = a 2 B a field and using the conformal-time variable (77 with 
r) = 1/a > and ' = d/dr/), the Fourier decomposition of ([3]) leads to 

B'( n) +n 2 B {n) = -2KB (n) , (4) 

since R a t = (2K/a 2 )h a i, to zero order. In the above expression, which does not explicitly depend 
on the matter component, K = 0, ±1 is 3-curvature index of the background Friedmann model 
and B a = Br n \Q^ is the harmonically decomposed magnetic vector, with D a B n = = Qa = 

D a Qa^ and D 2 Q^ = — (n/a) 2 Qi n \ The Laplacian eigenvalue takes continuous values, with 
n > 0, when K = 0,-1 and discrete ones, with n 2 > 3, for K = +1. When K = 0, the 
right-hand side of (Jl|) vanishes and we recover Eq. ([I]). Otherwise, we need to account for the 
effects of the background geometry. 

The curvature term on the right-hand side of expression ([!]) can in principle modify the 
adiabatic decay-law of the 5-field. For a K = — 1 background, in particular, 

/3' ( ' n) + (n 2 -2)£ (n) =0, (5) 

with n 2 > 0. On large enough scales, with n 2 < 2, the solutions to (5) no longer have the standard 
wave-like nature associated with the flat and the closed FRW hosts. These wavelengths include 
what one may regard as the largest subcurvature modes (i.e. those with 1 < n 2 < 2) and the 
supercurvature scales (having < n 2 < 1). Eigenvalues with n 2 = 1 correspond to the curvature 
length with physical wavelength A = a. Well inside this scale, the 3-space is practically flat, but 
beyond it the curvature dominates. Note that, although they are often omitted, supercurvature 
modes are necessary if we want perturbations with correlation lengths bigger than the curvature 
radius (see |10| for further discussion). Here we will focus on the largest subcurvature modes. 
Let us introduce, for convenience, the scale-parameter k 2 = 2 — n 2 with < k 2 < 2. Then, the 
largest subcurvature scales correspond to the range < k 2 < 1, while the interval 1 < k 2 < 2 
contains the supercurvature lengths. In the new notation, the solution of Eq. ([5]) reads 

B (k) = a' 2 [d cosh(|k|r/) + C 2 sinh(|k|r/)] = a~ 2 [C 3 e |k|j? + C 4 e~^} , (6) 

where the C^s are the integration constants. As we will show next, magnetic fields obeying the 
above evolution law can experience superadiabatic amplification, without modifying conventional 
electromagnetism and despite the conformal flatness of their FRW host. 

Suppose that the background model is a Milne-type universe: a vacuum, spatially open 
FRW spacetime with a = t. The latter immediately translates to e v oc a, which substituted into 
solution ([6]) leads to 

B{k) = C ia M- 2 + C 2 a-^- 2 , (7) 

with C\ 2 = constant. Consequently, all magnetic fields spanning lengths with < k 2 < 2 are 
superadiabatically amplified and their amplification strengthens with increasing scale. Close 
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Figure 1: The ratio d(ln5)/d(lna) for the dominant magnetic mode (vertical axis) versus the 
scale parameter k (horizontal axis), according to solution (|7|). Superadiabatic amplification 
occurs on all scales with < k < y2. The k = 1 value gives the curvature scale where B oc a -1 . 
Stronger amplification occurs on supercurvature lengths, with 1 < k < y/2. Thus, as k — » y/2 
and we approach infinite wavelengths, B oc a^~ 2 . At the k = limit the 5-field is well inside 
the curvature radius and the adiabatic decay is restored. 



to the curvature-scale threshold, that is for k — * 1, the dominant mode is B^ oc a , a rate 
considerably slower than the adiabatic a _2 -law. The latter is only restored in the k = limit, 
namely on small scales where the curvature effects are no longer important. Even stronger 
amplification is achieved on supercurvature scales, with Bqa oc a^ -2 at the homogeneous limit 
(i.e. as k — ► y/2 - see Fig. [I] for a summary). The Milne universe is probably the simplest, but 
not the only FRW background, that supports magnetic amplification of the superadiabatic-type. 
Solution ([6]) leads to similar results in open universes with p = p/3 as well. To verify this recall 
that the scale factor of a radiation-dominated Friedmannian spacetime with K = — 1 evolves 
according to a oc sinhr/ (e.g. see [H]). Then, focusing on the curvature scale for simplicity, 
Eq. © ensures that a magnetic field with k = 1 never decays faster than Bh\ oc a -1 and is 
therefore superadiabatically amplifiedH 

The Milne and the radiation-dominated Friedmann universe serve as very straightforward 
mathematical counterexamples, showing that the conformal flatness of the FRW spacetimes 
alone cannot guarantee an adiabatic decay for cosmological magnetic fields. For the purposes 
of physical cosmology, the Milne universe is the future attractor of (conventional) open FRW 
models. Thus, solution ([7J) could describe the late-time evolution of cosmological -B-fields. The 
radiation example, on the other hand, refers to the early universe and can drastically increase the 
expected strength of large-scale magnetic fields. The high conductivity of the pre-recombination 
plasma, however, means that we cannot a priori ignore the electric currents in Eq. ([3]). At 
the infinite conductivity limit, in particular, these currents are expected to dominate, eliminate 
the electric fields, and cause the magnetic component to freeze-in and dissipate adiabatically. 

3 One can easily verify that, on the curvature scale (i.e. for n 2 = 1), a magnetic field evolving like B a oc a -1 
is a solution of Eq. <(3j when p = and also when p = p/3. The former equation of state corresponds to a Milne 
universe and the latter to a radiation-dominated FRW model (with K = — 1 in our case). 
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Nevertheless, causality should confine the effects of the post-inflationary electric currents inside 
the horizon and therefore leave Eq. ([3]) unaffected on super-Hubble scales. Recall that our 
electromagnetic field crossed the horizon before the onset of the radiation era, during the poorly 
conductive epoch of inflation. 

To estimate the implications of the above described effects for primordial magnetic fields, 
it is important to note that large-scale -B-fields experience an analogous superadiabatic type of 
amplification in open FRW universes with an inflationary (i.e. p = —p) equation of state |12j . 
In these models, fields with coherence lengths close - and beyond - the curvature scale also 
experience superadiabatic amplification triggered by the same 3-curvature effects we described 
above. To be precise, -B-fields near the curvature scale were found to obey the evolution law 

B (1) = C 3 (1 - e 2 ") a~ l + C^a- 2 , (8) 

where rj < and 63,4 constants (see [H] for details). Thus, for most of the inflationary phase 
(i.e. as long as rj <C 0), we have Bn\ oc a" 1 and the field is superadiabatically amplified. The 
adiabatic decay rate is recovered only at the end of inflation, as rj — » 0" . We emphasise that the 
reduction in the magnetic decay-rate is possible because inflation in curved Friedmann models 
does not lead to a globally flat de Sitter space. Although inflation can dramatically increase the 
curvature radius of the universe, it does not change its spatial geometry. Unless the universe 
was perfectly flat from the beginning, there will always be a scale where the 3-curvature effects 
are importantly It is on these lengths that the -B-fields can be superadiabatically amplified. 

Following j!2j . a magnetic field that survived a de Sitter-type inflation in a FRW cosmology 
with K = -1 has energy density p B ~ HT 51 (M/10 17 GeV) 8/3 (TRtf/10 9 GeV)~ 2/3 A M 2 pc/ o 7 , and 
a current typical (comoving) strength from ~ 10~ 35 to ~ 10~ 33 Gauss, depending on the pa- 
rameters of the defining inflationary model. Note that M is the energy scale of inflation, Trh 
is the reheating temperature, p 7 is the radiation energy density and \Mpc is the current scale of 
the amplified -B-fieldH The latter is close to the curvature scale which, for a marginally open 
universe - with 1 — ft ~ 10 -2 today, lies between 10 4 and 10 5 Mpc. These scales are far larger 
than the minimum magnetic length required for the galactic dynamo to worklfl Nevertheless, 
once galaxy formation starts, the fluid motion should force the magnetic force lines to break 
up and reconnect on lengths similar to the size of a collapsing protogalaxy. Magnetic fields in 
the range of 10 -35 to 10~ 33 Gauss are far stronger than any other conventional -B-field that 
went through an epoch of inflation. So far, similar strengths have only been achieved outside 
standard electromagnetism. Moreover, seeds around 10 -34 G (or less) can sustain the galactic 
dynamo if our Universe is currently dark-energy dominated [13J. 

The above quoted strengths assume that the magnetic component freezes-in after the end 
of inflation and that the ratio Pb/p^ remains constant throughout the radiation epoch and 
until today [12]. Hence, the residual -B-field will increase further if it is superadiabatically 
amplified during the radiation era as well. Following our earlier discussion, this is possible 

4 In probabilistic terms, an exactly flat universe is a set of measure zero and therefore physically unrealisable. 

5 Assuming M ~ 10 17 GeV and Ho ~ 70 km/sec Mpc, one obtains a residual B-field of the order 10 -35 G when 
Trh ~ 10 9 GeV. A reheating temperature close to 10 3 GeV, on the other hand, leads to B ~ 10~ 33 G. In general, 
the higher the scale of inflation and the lower the reheating temperature the stronger the amplification [12] . 

6 The dynamo mechanism requires magnetic seeds with a (collapsed) coherence length close to 100 pc. This 
corresponds to a comoving (pre-collapse) scale of approximately 10 Kpc. 
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for superhorizon-sized magnetic fields of inflationary origin, because they are not affected by 
the electric currents of the post-inflation universe. In such evolution law of a 

during the whole of the radiation epoch, will add several orders of magnitude to the residual 
-B-field. For example, the above quoted magnetic seed of 10~ 35 G will 'grow' further during 
the radiation era. This field, which has length close to the curvature radius and corresponds to 
Trh ~ 10 9 GeV, will be superadiabatically amplified by TRH/T rec ~ 10 19 orders of magnitude 
and reach a comoving strength of up to 10~ 16 G. Note that, when 1 — Qq ~ 10 -2 , magnetic 
fields spanning lengths near the curvature scale of the universe, remain outside the horizon 
(and therefore are superadiabatically amplified) throughout the radiation era. Also, although 
they leave primordial nucleosynthesis and the Cosmic Microwave Background unaffected |14j . 
magnetic seeds of the order of 10 -16 G are astrophysically important because they can sustain 
the galactic dynamo even in conventional universes with zero dark energy [15]. With this in 
mind, a more detailed (most likely numerical) analysis is necessary to establish the full spectrum 
of the residual B-field. Here we have confined ourselves to effects close to the curvature scale, 
which probably means that the estimated magnetic strengths are the maximum possible. 

Large-scale magnetic fields have been observed almost everywhere in the universe. From 
galaxies and galaxy clusters, to super-clusters and remote high-redshift protogalactic structures, 
observations have repeatedly detected coherent S-fields of micro-Gauss strength. Despite this 
widespread presence, however, the origin of cosmic magnetism is still unknown and a matter of 
debate. Inflation has long been seen as our best mechanism for generating large-scale magnetic 
fields, because it naturally leads to superhorizon-sized correlations. The resulting fields, however, 
were always considered too weak to have any astrophysical significance. This was attributed 
to the conformal invariance of Maxwell's theory and to the conformal flatness of the FRW 
spacetimes. The two were believed to guarantee an adiabatic decay for any cosmological B- 
field, irrespective of the electric properties of the universe. 

Contrary to this widespread belief, however, magnetic fields on FRW backgrounds do not 
always experience an adiabatic, Minkowski-like depletion. Departures from the a _2 -law, with- 
out breaking away from conventional electromagnetism, are possible in FRW models with 
non-Euclidean spatial hypersurfaces, because those spacetimes are not globally conformal to 
Minkowski space. This geometrical subtlety, which has long been known within the relativity 
community, but has been largely overlooked in studies of cosmological magnetic fields, is central 
to our report. Together with the vector nature of the Maxwell field, which brings the 3-curvature 
into play, it is essential in understanding the behaviour of large-scale magnetic fields in FRW 
cosmologies and avoids the need to deviate from Maxwellian electromagnetism in our quest for 
astrophysically relevant cosmological magnetic seeds. 
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